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ABSTRACT 

We present an analysis of optical and ultraviolet Hubble Space Telescope photometry for evolved 
stars in the core of the distant massive globular cluster NGC 2419. We characterize the horizontal 
branch (HB) population in detail including corrections for incompleteness on the long blue tail. The 
majority of the horizontal branch stars can be identified with two main groups (one slightly bluer 
than the instability strip, and the other at the extreme end of the HB). We present a method for 
removing (to first order) lifetime effects from the distribution of HB stars to facilitate more accurate 
measurements of helium abundance for clusters with blue HBs and to clarify the distribution of stars 
reaching the zero-age HB. The population ratio R — Nhb/Nhgb implies there may be slight helium 
enrichment among the EHB stars in the cluster, but that it is likely to be small {AY < 0.05). An 
examination of the upper main sequence does not reveal any sign of multiple populations indicative 
of helium enrichment. 

The stellar distribution allows us to follow how the two main types of stars evolve after the HB. 
We find that the transition from stars that reach the asymptotic giant branch to stars that remain 
at high temperatures probably occurs among the extreme horizontal branch stars (EHB) at a larger 
temperature than predicted by canonical evolution models, but qualitatively consistent with helium- 
enriched models. Through comparisons of optical CMDs, we present evidence that the EHB clump 
in NGC 2419 contains the end of the canonical horizontal branch, and that the boundary between 
the normal HB stars and blue hook stars shows up as a change in the density of stars in the CMD. 
This corresponds to a spectroscopically-verified gap in NGC 2808 an d an "edge" in cu Cen. The more 
clearly visible HB gap at ~ 23.5 identified bv lRipepi et al.l (|2007D appears to be too bright. 

Once corrected for lifetime effects, we find that NGC 2419 is currently converting about 25 - 31% 
of the first-ascent red giant stars in its core into extreme blue horizontal branch stars — the largest 
fraction for any known globular cluster. A comparison of upper red giant branch with theoretical 
models indicates there is a slight deficiency of bright red giant stars. This deficiency occurs far 
enough below the tip of the red giant branch that it is unlikely to be associated with the production 
of extreme horizontal branch stars via strong mass loss before the core helium flash. 
Subject headings: stars: horizontal branch — stars: evolution — stars: luminosity function — stars: 
mass loss — globular clusters: individual (NGC 2419) 



1. INTRODUCTION 

NGC 2419 is the fourth brightest gl obular cluste r 
known in the Milky Way (My = -9.58; lHarrislll996l ). 
NGC 2419's large mass would normally make it ideal for 
studying short phases in the lives of evolved low-mass 
stars because massive cluster (having more stars) should 
produce a greater chance of catching some stars in those 
phases. Until recently though, it has been neglected be- 
cause it resides far in the outer halo of the Milky Way 
(approximately 90 kpc from the center) . Massive clusters 
can make it possible to observe stars in short-lived evolu- 
tion phases because there is a greater chance of catching 
some star in such a phase when there are more chances. 

Evolved stars have unusually strong influences on the 
luminosity-integrated properties of a galaxy: in old stel- 
lar populations, the overall color and the ultraviolet emis- 
sion result from them. The luminous giants at the tip of 
the red giant branch (TRGB) are frequently used as stan- 
dard candles in resolved stellar populations. Extremely 

Electronic address: erics@sciences.sdsu.edu 

^ Based on observations with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astron- 
omy, Inc., under NASA Contract NAS 5-26555. 



blue HB (hereafter, EHB) stars are a leading candidate 
for galactic ultraviolet emission. These populations are 
linked at least in an evolutionary sense, but it remains a 
longstanding problem to explain the details of how HB 
stars (and more specifically the EHB stars) are produced. 
There has been much recent interest in multiple stellar 
populations having different chemical compositions (par- 
ticularly helium, but also heavier elements) within indi- 
vidual clusters (N GC 280 8.lD'Antona et a l. 2005: uj Cen, 
iBedin et al][200l M13. ICaloi fc D'Antoia 2005: NGC 



6218. lCarretta et al.ll2007t NGC 6388. lBusso et al.l l2007l: 
NGC 6441. ICaloi fc D'A^tom] [20071) However, even al- 
lowing for composition effects, it still seems to be nec- 
essary for giant stars to lose most of the mass in their 
envelopes if EHB stars are to be produced — the turnoff 
mass of a helium-enriched population cannot be reduced 
enough to have a giant star consume its envelope before 
it reaches the TRGB. 

In this paper, we examine deep photometry of the 
evolved stars in the core of the cluster using high res- 
olution Hubble Space Telescope (HST) imagery. We 
have three main goals: to characterize brief evolutionary 
phases and their effects on the integrated light of the clus- 
ter, to examine the paths (often brief) that lead from one 
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evolutionary phase to another, and to use the stellar pop- 
ulations to constrain chemical evolution within the clus- 
ter. EHB and bright red giant branch (RGB) stars are 
particularly important to the integrated colors of clus- 
ters. EHB st ars can be seen in several earlier CMDs 
for NGC 2419 (lHarrislll997l : lStetsonll200^ . although the 
impression was that the popul ation was small due to in - 
completeness at the faint end. iDalessandro et all ()2008D 
recently presented a reduction of a portion of the avail- 
able archived HST data, for the cluster that clearly shows 
the faint end of the blue HE tail, and that the EHB pop- 
ulation is a substantial fraction of all HB stars in the 
cluster. In section 13.11 we conduct a deeper analysis of 
the HB population of the cluster. 

The details of the transition from the RGB to the HB 
remain poorly und erstood. In the case of the massive 
cluster NGC 2808, iSandauist fc Martell (|2007D found a 
deficit of bright red giants compared to theoretical pre- 
dictions that may be linked to the EHB stars there — if 
a star loses enough mass on the red giant, it can leave the 
red giant branch before the helium flash at the TRGB, 
igniting helium later in a "hot flash" that ultimately de- 
posits it at the blue end of the HB. NGC 2419 appears to 
be more efficient at producing EHB stars than most clus- 
ters, so in ij3.3l we examine the upper RGB for signatures 
of EHB star production. 

Integrated luminosity has been linked to a number 
of unusual groups of stars within globular clusters, in- 
cluding blue stragglers (their relative frequency is anti- 
correlated with total luminosity), extremely blue hori- 
zontal branch stars, and populations with different age 
and/or chemical composition. To date, NGC 2419 has 
shown little indication that it contains anything but 
a stellar population with a single composition ([Fe/H] 
= -2.1; ISuntzeff et al|[l988! ) and age. However, NGC 
2419 is now known to have a strongly double-peaked HB, 
and so we examine the helium abundance indicator R in 



2. OBSERVATIONS AND DATA REDUCTION 

This study has focused on the evolved stellar popula- 
tions of the cluster core. The archival HST datasets that 
were reduced are listed in Table [U and a map of the fields 
is shown in Fig. [TJ Images from the ACS WFC and HRC 
instruments were processed using the DOLPHOT pho- 
tometry package^ with its module tuned for ACS data. 
Individual frames were obtained from the HST archive in 
order to get photometry on the giant stars. The WFPC2 
images were analy zed using the HSTPhot^ photometry 
package ()Dolphinl[200d ). Although different filters were 
used for imaging in different epochs, the fields overlapped 
in many cases, and in each field we have photometry in 
the F555W and F814W filters. 

An RA-DEC offset coordinate system was set up us- 
ing the WFPC2 datasets. Positions a.re giv en relative 
to the cluster center given by ,Harri9 (|1996f ). This po- 
sition agrees extremely well with the dete r minat ion of 
the center of gravity bv IDalessandro et al.l ()2008( ). Im- 
ages in the ACS system are subject to a lot of geometric 
distortion due to the off- axis positions of the cameras, 
so we used fourth-order geometric transformations ob- 

^ http;//purcell. as. arizona.edu/ andy/dolphot/ 
^ http:/ /purcell. as. arizona.edu/hstphot/ 
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Fig. 1. — Sky positions of HST fields that were analyzed in 
this p aper. Positions are given in arcseconds relative to the Harris 
(119961) position given for the cluster center. WFPC2 fields are 
identified with the proposal ID number (see Table [l} except where 
the WFPC2 field almost completely overlaps another WFPC2 field. 
The deepest photometry was taken in the fields having bold lines. 
The ACS HRC field is the small rectangle surrounding (0,0). 

tained from the most recent Image Distortion Correc- 
tion Table (IDCTAB) file (qbul641sjJdc.fits)4 in order 
to convert pixel positions into relative sky positions. Co- 
ordinates for stars in the WFPC2 images were derived 
from the METRIC task in the IRAF package STSDAS. 
METRIC positions typically have absolute position er- 
rors that are a significant fraction of an arcsecond, but 
we are primarily interested in relative positions. Af- 
ter initial identifications of stars were made in the dif- 
ferent datasets, preliminary coordinate transformations 
were derived for all of the WFPC2 and ACS fields. The 
coordinate transformations were then refined with addi- 
tional cross-identifications, and these were used to match 
up all stars measured in the datasets. 

The HST datasets were analyzed to maximize the de- 
tection of different populations of evolved stars in the 
cluster. Identification of horizontal branch stars is gen- 
erally more effective in shorter wavelength filters when 
most of the HB stars are bluer than the RR Lyrae insta- 
bility strip, as in NGC 2419. In this case, however, the 
optical CMDs from the ACS WFC field were deep and 
had high spatial resolution. At the same time, the op- 
tical CMDs were significantly incomplete in the core of 
the cluster even with the high resolution. In Fig. [21 we 
show the optical WFC CMD for samples in and outside 
60" from the cluster center. From artificial star tests (see 
t j2.ip , we find that the WFC data is more than 80% com- 
plete even down to the faint end of the EHB for r > 60". 
Because they also covered the widest area, they became 

* http: / / www.stsci.cdu/hst / obscrvatory/cdbs/SIfilcInfo / ACS/ ACSDistortion' 
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TABLE 1 
HST Datasets for NGC 2419 



Prop. ID'' 


P.I. 


Instrument 


Filters 


5481 


Hesser 


WFPC2 


F555W, F814W 


7628 


Whitmore 


WPPC2 


F336W, P380W, P439W, F450W, F555W, F675VV, F814W 


8095 


Ibata 


WFPC2 


F555W 


9601 


Koekemocr 


WFPC2 


F300W, F336W, F439W, F450W, F675W 


7630 


Casertano 


WFPC2 


F300W, F555W, F814W 


6937 


Stiavelli 


WFPC2 


F555W 


8821 


Hiess 


WFP(J2 


F555W, F814W 


9043 


Tonry 


WFPC2 


F555W, F814W 


9591 


Whitmore 


WFPC2 


F814W 


9666 


Gilliland 


ACS WFC 


F435W, F555W, F814W 


10815 


Brown 


ACS HRC 


F250W 



^ Observations with almost completely overlapping fields are grouped together. 



our primary choice for selection. Had the ultraviolet im- 
ages been deeper, they would have been more reliable for 
detecting the bluest HB stars. The ACS HRC frames in 
the F250W filter covered a small field centered about 12" 
from the cluster core, allowing us to confidently identify 
11 faint HB stars that would otherwise have been rejected 
by quality cuts that were applied to the next highest res- 
olution datasct (the ACS WFC images). Fig. [3] shows 
CMDs of HB stars detected in both the HRC and WFC 
images. Based on these considerations, we believe the HB 
star sample in the HRC field is virtually 100% complete. 
Observations in the F336W filter of the proposal 7628 
field also provide good detections of bright EHB stars, 
although the detection limit brighter than the faint end 
of the EHB (see Fig. [H]). The situation in the F300W 
filter was similar, so they were only used for verification 
of brighter candidates. 

Because some of the datasets have image sequences 
that are suitable for searching for variables, we checked 
the RR Lyrae popul ation in the H ST fields. We com- 
puted the J statistic (|Stetsonlll996l ) for the combination 
of F555W and F814W images from the proposal ID 6937 
and 7630 field, and found that this efficiently identified 
variable stars in the field. We confirmed most of the 
iPinto fc Rosinol (|1977D variables that fall in the field (V2, 
V3, V4, V6, Vll, V13, V26, V29, V30, V33, V35, V38, 
V39, V40) and have detected others. The coordinates 
given for V37 were inconsistent with the ' Pinto fc Rosinol 
finding chart position, which puts it outside the HST 
fields. The coordinates for V38 and V40 were both sig- 
nificantly in error, but the stars were identified using the 
finding chart and their variability was confirmed. We 
also were unable to find an obvious counterpart to their 
V41. Most of the RR Lyrae stars could be tentatively 
identified in the WFC datasets using their CMD posi- 
tions, but there were 4 stars that ended up in the middle 
of the RGB as a result of the sampling of the light curve 
with a small number of observations. 

The CMDs used for the selection of RGB stars are 
shown in Fig. |4l In the RGB star sample, the primary 
contaminant is AGB stars. For this reason we also re- 
lied heavily on CMDs with short wavelength filters in 
order to separate AGB and RGB stars by temperature. 
CMDs with ultraviolet observations on the magnitude 
axis can help turn temperature (color) differences into 
luminosity (magnitude) differences, which are often eas- 
ier to detect. Our primary choice for RGB identification 
used the F336W filter from the proposal 7628 field be- 



cause this most cleanly and reliably separated AGB stars. 
As our second choice we used ACS WFC observations in 
the F435W filter and WFPC2 observations in the F439W 
and F450W filters. Stars in the longest WFC exposures 
were unfortunately saturated, although we were able to 
use saturated star photometry from DOLPHOT in most 
cases. (We have not included saturated star photometry 
in the averaged photometry described later.) If there was 
some ambiguity in identifying a star as RGB or AGB, it 
usually came in the CMDs due to (relatively) small color 
separation between RGB and AGB. The AGB clump can 
easily be seen in 19.9 < Fi35W < 20.5, but the AGB se- 
quence mostly merges with the RGB brighter than that. 
Photometry in the F300W (WFPC2) and F250W (HRC) 
filters were also useful, but again, we found that they 
were sometimes unreliable toward the tip of the giant 
branch due to measurement uncertainties. Toward the 
faint end of the sample though, there is a large and use- 
ful magnitude difference between the RGB and the AGB 
clump. 

2.1. Artificial Star Tests 

In order to fully characterize the faint HB population, 
we conducted a series of artificial star tests on the ACS 
WFC frames, which comprises the great majority of the 
surveyed field. Before we ran the tests, we first defined 
a fiducial line running through the blue HB tail using 
the F435W filter for the magnitude axis, and using both 
i^435W^ - F555W and i^435W^ - F8UW colors. To en- 
sure that the two artificial star colors were consistent, 
we used the photometry for a select group of stars to de- 
fine the fiducial. Artificial star magnitudes were selected 
randomly with a flat luminosity function, and colors were 
generated from quadratic interpolation between the fidu- 
cial points. 

In each of the artificial HB star tests, stars were placed 
in a grid of cells 30 pixels wide. The stars were ran- 
domly placed within cells, and the grid itself was shifted 
from test to test. A total of 5 runs containing 91926 
stars were completed. Stars were identified as recov- 
ered if they were detected in the F435W and F814W 
images and passed cuts on sharpness (|SHARP| < 0.5), 
(< 5), and crowding. With regard to crowding, the 
photometry for a star is often affected by light from 
nearby neighbors. Photometry was derived from pro- 
file fitting in which overlapping stellar images are fit si- 
multaneously with point-spread functions. The crowding 
measure involves the resulting difference in brightness if 
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Fig. 2. — ACS WFC color-magnitude diagrams with horizontal branch stars passing criteria for detection in our artificial star trials (blue 
crosses; see i|2.1l l. divided by projected distance from the cluster center. 



nearby neighbors were not simultaneously fitted. Stars 
were rejected if they would have been more than 0.3 mag 
brighter if nearby neighbors had been neglected. (While 
this probably rejected some real stars, if we apply the 
cuts consistently to the artificial stars and the real stars, 
we can derive accurate incompleteness corrections.) 

For the analysis, stars were divided into annuli (20" 
wide; approximately one core radius) and magnitude (0.2 
mag wide) bins. According to the artificial star tests (see 
Fig. [5]), the bright HB sample is almost 100% complete 
throughout the HST fields. The worst incompleteness 
occurs in the cluster core, where EHB stars in the WFC 
field are only about 29% complete within one core radius. 
However, we used UV observations to improve the recov- 
ery of EHB stars in the densest part of the core. Because 
crowding is negligible in the HRC images and the faint 
end of the HB is well above the faint limit of the images, 
we have assumed that our HB sample is 100% complete 
in the HRC field. 

2.2. Photometry 

With such a heterogeneous set of images and filters, 
we spent some effort to merge the photometric datasets 
together where possible. We emphasize that the merged 
photometry was not used in categorizing stars by evolu- 
tionary state — that was done using individual datasets 
because the relative precision was higher. We used the 



observations of the WFPC2 proposal 7630 field as the ref- 
erence because it showed the greatest amount of overlap 
with the other fields. (We found significant zeropoint 
differences between the two ACS WFC chips, which is 
why they were not used.) The photometry that could be 
merged in this way used the F300W, F555W, and F814W 
filters. 

In the F555W filter, we found no signs of color depen- 
dence among the observations in the different datasets, 
so we only corrected for zeropoint differences. The off- 
sets (labeled Cq) are provided in Table [U and are in 
the sense of (7630 - dataset). Residuals versus F555W 
magnitude are plotted in Fig. [6l In the F814W filter, 
we found a significant color dependence in the residu- 
als between the 5481 and 7630 photometry. Note that 
all of the WFPC2 photometry was rereduced using the 
same software and photometry parameters, so this differ- 
ence appears to come from instrumental effects. Color- 
dependent terms are labeled ci in the table. Residuals 
versus F814W magnitude are plotted in Fig. [7l and the 
dependence on color is shown in Fig. [51 We note that the 
color dependence of the transformation from the WFC 
to the WFPC2 in t he F555W and F814W filters evalu- 
ated bv lSirianni et a l. (2005) are consistent with zero (to 
within 1 a) for all color choices. 

We found a fairly large zeropoint difference and color 
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TABLE 2 

Photometric Transformation to Proposal 7630 
Photometry 



Fig. 3. — Color-magnitude diagrams showing the selection of extreme HB stars (blue crosses) from the ACS HRC F250W filter, and the 
corresponding stars in the ACS WFC field. Triangles sh ow p robable evolved HB stars in the ACS HRC field, with the triangles pointing 
roughly in the direction of their expected evolution (sec i|3.2l l. 

zeropoint differences in comparison to tfie flight photo- 
metric standard system.) The difference in F300W mag- 
nitudes as a function of magnitude and color are shown 
in Fig. El 

The WFC F435W and WFPC2 F439W filters are also 
close in wavelength coverage, although the NGC 2419 
field was only observed in these filters using WFPC2 in 
proposal 7628. We compared the photometry, and found 
a small color term is able to do a good job of transforming 
the two. The color term (cq = 0. 030) is slightly large r 
than the one given in Table 20 of lSirianni et al.l (120051 
but does a better job at simultaneously bringing the HB 
and RGB populations into consistency. Fig. [10] shows 
the comparison between the WFPC2 and transformed 
WFC data. 

We have generally left our photometry in the VEGA- 
MAG flight system. When transformations to UBVRI 
bandpasses are needed for comparisons to external pho- 
tometry for ot her clusters, we have used transformation 
equations from lSirianni et all ()2005D for ACS WFC im- 
ages, and the most recent transformations maintained 
by A. Dolphin^ for WFPC2 images. The flight system 
photometry is presented in Tabled] 

3. ANALYSIS 
^ |http: //purcell ■ as . arlzona. edu/wf pc2_calib/| 



Prop. ID 


F300W F555W 

Co Cl Co 


F814W 

Co Cl 


5481 


-0.068 


-0.141 0.0778 


7628 


0.113 0.0201 -0.017 


-0.019 


8821 


-0.005 


0.004 


9591 (S)'^ 




-0.005 


9591 (N)^ 




-0.032 


10815 (WFCl-E)'' 


-0.033 


0.011 


10815 {WFC2-W)'^ 


-0.078 


-0.037 



The two proposal ID 9591 fields are labeled south or north based 
on their positions, and both are shown in Fig. [T]'' The fields of the 
two ACS WFC chip fields in proposal ID 10815 fields are labeled 
east or west based on their positions, and are shown in Fig. [T] 

dependence of the photometric residuals in the F300W 
filter. Time-depende nt zeropoint changes a. re well-known 
in ultraviolet filters (jHoltzman et al.lll995[ ). but the be- 
havior in F300W has not been thoroughly mapped out 
because the filter is infrequently used. We have therefore 
only corrected for relative differences between the two 
F300W datasets, again using the proposal 7630 photom- 
etry as the standard. (Of all of the filters considered here, 
the ultraviolet filters are most likely to have significant 
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Fig. 4. — The color-magnitude diagrams used to select red giant branch branch (red) and asymptotic giant branch (green) stars in the 
HST fields. No cuts on position with the cluster were applied to any of the CMDs. 
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Fig. 5. — Completeness measures for artificial horizontal brancfi 
stars as a function of F814W magnitude and distance from the 
cluster center. The lines are identified with the radius of the center 
of the annulus (all annuli are 20" wide). 

3.1. The Horizontal Branch 

The HB morphology has been recently described 
(based solely on ACS photometry) in iDalessandro et al.l 
(|2008i), but the focus of that study was on the blue 
stragglers. Here we attempt to fully characterize the 
evolved populations, and in that vein we have identi- 
fied 994 HB stars in the union of the HST fields. Our 
F435M^ and i^814W^ magnitudes from the ACS data 
are approx imately equivalent to the B and / magni- 
tudes from IDalessandro et all . For most of the discus- 
sion below, the two pairs can be used almost inter- 
changeably (color t ransformation coefficients are small; 
ISirianni et a'Lll2005f ). There are a number of interesting 
features of the CMDs deserving of comment. 

The HB has an extremely large total extent in effective 
temperature, although the red end of the HB and the 
RR Lyrae instability strip are sparsely populated, as is 
theoretically expected for the cluster's metallicity. There 
are at least four clear breaks in the distribution of stars 
that can be seen in Fig. [TT] 

1. B — I = 0.7: This break approximately corresponds 
to the blue end of the instability strip, but there is 
a clear drop off in the number of stars redder than 
this point. Although we have been able to iden- 
tify RR Lyrae variable stars from the datasets, we 
generally do not have good enough coverage of the 
light curve to derive precise average magnitudes. 

2. B « 21): while not a gap in the distribution, this 
marks the blue edge of the most heavily populated 
part of the HB. The HB is more sparsely populated 
out to the next fainter break. 

3. B w / « 22.25: This is the brighter of the two 
clear gaps in the HB distribution. We will refer to 
the stars with 22.25 < B < 23.5 as the inter-gap 
population. 



4. S « 23.5: The fainter of the two gaps marks 
the bright edge of the collection of extreme HB 
(EHB) and blue hook (BHk) stars identified by 
IDalessandro et all . Those authors tentatively iden- 
tified this as a gap. 

IDalessandro et all compared the HB distribution with 
the more massive cluster uj Cen, and found similarities 
in the overall distribution — notably comparable per- 
centages of stars in the EHB/BHk group. The compar- 
ison with {B, B ~ I) data for M13 in Fig. [11] clearly 
shows that NGC 2419's HB extends approximately 1 mag 
fainter than the greater part of the EHB clump for M13 
(although a handful of stars have been identified fainter 
than the clump). A comparison with two other clusters 
with long blue HB tails (NGC 2808 and M15) is shown 
in Fig. [121 In each case, the number of HB stars appears 
to sharply drop at approximately the same magnitude 
(relative to the edge of the RR Lyrae instability strip): 
B ~ 21 in NGC 2808, and B - 25.1 in NGC 2419. M15 
is a cluster of similar [Fe/H] and comparable but lower 
mass that converts only a has a tiny proportion of its 
stars into EHB stars. NGC 2808 has a mass that is con- 
siderably closer to that of NGC 2419, but it has a much 
higher [Fe/H]. It produces more EHB stars than M15, 
but the EHB stars are still a small fraction of the total 
HB stars. 

The majority of the cluster's HB stars fall in a compact 
group in the CMD (20.6 < B < 21, 0.2 < B - / < 0.6; 
see Figure ITT]) , with less dense extensions to the red and 
blue. There is a fair amount of evidence to suggest that 
all or nearly all of the stars redward of the main peak 
have evolved to their present colors from the main peak. 
The presence of AGB stars (especially ones in a "clump" 
at the base of the AGB) implies redward evolution af- 
ter the HB phase, and at least some models (see the 
next paragraph) predict stars should spend a substantial 
time redward of their starting point on the HB. Proba- 
bly the strongest evidence comes from the comparison of 
the period-amplitude diagr am for NGC 2 419 RR Lyrae 
stars with that of M3 (Ripc pi et aTll2007f ). M3 is a clus- 
ter which has strong populations of stars both redder 
and bluer than the instability strip, and so it is legiti- 
mate to expect from theoretical evolutionary tracks that 
there are relatively unevolved RR Lyrae stars as well as 
stars that have evolved into the instability strip from the 
blue HB. The lower average densities of the evolved stars 
are expected to result in larger periods, and this is sup- 
ported by the appea rance of two sequenc es in the period- 
amplitude diagram (ICac ci ari et"alll2005f) . NGC 2419 RR 
Lyraes fall nearly exclusively am ong the "evolved" p op- 
ulation of M3 stars (see Fig. 3 of lRipepi et al.ll2007t l. In 
this context, star V2 appears to be an unusually bright 
evolved RR Lyrae star. The CMD position identifies it as 
a "supra-HB" star (see §3.2p . We have identified V2 from 
the finding chart published by Pinto & Rosino (1977), 
and from our own detection of its variability in the HST 
data. Our photometry implies that it was observed near 
minimum light in the F814 W filter (see Fig. [T3| . The 
only published light curve (jPinto fc Rosinolll977i ) implies 
it has a large amplitude (about 1.3 B magnitudes) and 
period, although their photometry indicates that it is 
comparable to other known RR Lyrae stars oscillating 
in the ftmdamental mode. At the same time, their pho- 
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tometry implies that it is also similar in brightness to 
other RR Lyrae s as well. This contradiction between the 
iPinto fc Rosinol photometry and ours suggests that this 
star should be studied in more dep th. 

HB sequences (e.g. Dotter et al. 2007f) imply that zero- 
age horizontal branch (ZAHB) stars for NGC 2419's low 
metallicity and primordial helium abundance always fall 
bluer than the instability strip, even for stars with no 
mass loss. Synthetic HB models also predict a signifi- 
cant number of variable and red HB stars for a zero-age 
population that begins at the position of the main peak 
(see Fig. [20|) . Blueward evolution in H B stars tends 
to oc cur in stars with lower initial Lhc/L (|Iben fc RoodI 
flOTOl) . where Lhc is the luminosity released by helium fu- 
sion. This means that the hydrogen-fusion shell rapidly 
adds mass to the helium core of the star. This can 
result from higher helium abundances or metallicities 
(both affecting the structure of the envelope, and pro- 
ducing higher shell fusion temperatures), but in general 
it requires a substantial hydrogen envelope in order to 
maintain the con ditions ne eded for a strong hydrogen- 
fusion shell (iSweigartlll987[ ) . When the hydrogen shell is 
weakened (by low envelope mass or low envelope helium 
abundance, for example), the color evolution is driven by 
the depletion of helium in the core as the star becomes 
more giant-like. So the evolution of HB stars can bear 
on questions of helium abundance, at least tangentially. 
The relatively low dispersion in magnitude among the 
red HB stars is consistent with the idea that the stars 
in the main peak have low (approximately primordial) 
helium abundance. An increase in Y of as little as 0.05 
produces an initial blue loop in the star's evolution, and 
the subsequent evolution toward the giant branch takes 



place at a significantly greater luminosity. (See Figure 
9 of iBusso et all l2007i for a comparison of much more 
metal-rich tracks with different helium abundances.) 

The more important issue is what can be learned about 
the population of the second peak on the EHB, as it 
is the bluer population of HB stars that is identified 
with helium-enriched gas when multiple populations are 
present. An understanding of this group is complicated 
by poorer sensitivity to temperature of the colors avail- 
able here, and the fact that some of these stars appears to 
exist beyond the end of the canonical HB, as seen in Fig. 
[TTl Because these stars are at or beyond the end of the 
HB, they are called extreme HB or blue hook stars (we 
will refer to them as the EHB/BHk group). There are 
major uncertainties in models of these stars: for example, 
it is still unclear exactly how the stars are formed and 
whether the formation mechanism has a strong effect on 
the position in the CMD, but the combination of flash 
nucleosy nthesis and convec tive mixing makes modeling 
difficult (H rown et aLll2001h . Regardless, in the analysis 
below, we first attempt to identify as many of the stars 
associated with this group as possible, and then examine 
the evolution paths they take. Because the evolution is 
seen most clearly in some of the ultraviolet datasets, we 
discuss it in ij3.2l 

In the ACS photometry (see Fig. [2]), we consistently 
see a tail of stars stretching from the EHB/BHk group 
toward the faint end of the blue straggler group and 
main sequence turnoff (23.6 < < 25, -0.3 < 

~ F8UW < 0.5). These are probably stars that 
are optical blends of an EHB/ BHk star with a main 
sequence star. Castell ani et al.l (2006) saw similar fea- 
tures in their optical CMD of NGC 2808. The features 
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posal 7628) and transformed ACS WFC F435W photometry in the 
sense of (WFPC2 - WFC). The solid line shows the median resid- 
ual. Zero point differences between the WFC chips were less than 
0.004 mag, and were not corrected. 

were not as strong there, probably because NGC 2808 
has a much smaller EHB population producing such ob- 
jects in the CMD. Fig. [Ml shows loci of blends of three 
EHB/BHk stars (marked by the faintest 3 x symbols) 
with main sequence an d giant branch stars [represented 
by a lDotter et al.ll2007l (hereafter, DSEP) isochrone hav- 
ing age of 13 Gyr and [Fe/H]= —2.1 that was fitted to 
the stellar locus]. Blends with turnoff stars fall near the 
blue straggler region near the kinks in the curves. While 
such blends could masquerade as blue stragglers in op- 
tical CMDs, they can be clearly identified in ultraviolet 
CMDs due to the much higher surface temperatures of 
the HB star. 

We use the results of the artificial star tests (t j2.ip to 
correct the distribution of detected HB stars for incom- 
pleteness. We took the total sample of identified HB stars 
and only kept stars that were detected in either the ACS 
WFC or HRC photometry and passed the same quality 
cuts as the artificial stars. This sample was made up of 
856 stars, and with incompleteness corrections, the total 
sample expected in the field was 954. The top panel of 
Fig. [T51 shows the detected and incompleteness-corrected 
HB distributions for the ACS fields. Based on this, the 
EHB population makes up more than 38% of the total 
HB population. 

3.1.1. Accounting for HB Star Lifetimes 

Theoretical models uniformly predict that stars at the 
blue end of the HB have 20-30% longer HB lifetimes than 
stars near the instability strip, so we need to correct for 
this effect to determine how frequently EHB stars are 
being produced compared to the rest of the HB popula- 
tion. Because the HB population is used in determining 
the initial helium abundance, corrections for differences 
in HB star lifetimes can make the measurements more 
accurate. Below we describe a method for correcting for 
variations in HB lifetimes. As much as possible, we wish 
to produce a correction that i) is independent of chemi- 



cal composition so that composition variations or imper- 
fectly known compositions do not introduce systematic 
errors, and ii) is based on observational quantities. 

For stars with the same initial compositions, we find 
that the lifetimes vary in similar ways as a function of 
stellar mass and effective temperature, although the ab- 
solute value of the lifetime does depend on composition. 
(In the following we have defined the end of the HB 
phase as the point when the central helium abundance 
Yc = 0.05; see the lefthand panels of Fig. [Tni Because 
HB stars evolve somewhat in temperature during this 
phase, we have used the temperature of the star when 
it is halfway through its HB lifetime as a representative 
Teff. We can imagine more complex algorithms for deter- 
mining this temperature, but even the use of the T^s at 
zero-age produced minor effects on the results.) Based 
on this behavior, we defined a weighting factor for each 
star in the HB sample: 

^ tga (log Te// ^3.85) 

tHB{\ogT,ff) 

The largest portion of the variation in HB lifetimes oc- 
curs on the blue tail for T^-ff > 10^ K in stars with low- 
mass hydrogen envelopes. I3y choosing logTe// — 3.85 
(near the blue edge of the instability strip) as the nor- 
malization point, we can correct the blue HB lifetimes 
to values representative of variable and red HB stars, 
which are common in globular clusters. The variation 
of the weighting factors with T^f / is shown in the right- 
hand panels of Fig. There is very little variation in 
the weighting factors with variations in [Fe/H], so that 
residual uncertainties are at the level of a few percent. 
Further, the weightings only remove the effects of posi- 
tion on the HB, and do not reference the absolute value 
of the HB lifetime, which does depend on composition 
and physics inputs to the stel lar evolution codes (e.g., 
the i2c(a, 7)1^0 reaction rate: iCassisi et al.ll2003D . 

The weighting factors also vary little with Y for 
log Toff < 4.2, although the agreement gets worse at 
higher Toff. Higher helium abundance produces HB stars 
whose lifetimes increase more rapidly with decreasing 
Toff than stars with lower helium. So for clusters with 
the majority of their HB stars having Toft ^ 16000 K, 
the weighting factors are almost entirely independent of 
composition, and corrected HB star totals can be used to 
gauge helium abundance directly using the R ratio. For 
clusters like NGC 2419 that instead have a larger frac- 
tion of their stars at higher Tefi, this is not as straight- 
forward, but the R value can still be used as a test of 
whether clusters stars have enhanced helium. In fact, if 
weighting factors having a canonical value of F « 0.25 
are used, then helium-enriched stars on the blue tail will 
be undercorrected for their long lifetimes. This has the 
effect of leaving the corrected number of HB stars and 
the R ratio too large, and implying a high Y. 

The weighting has the most utility when dealing with 
stellar populations that have the same composition. Al- 
though NGC 2419 has not shown any evidence of large 
star-to-star composition ([Fe/H], Y) differences in its 
CMD (via large dispersion in its evolutionary sequences 
or multiple populations), other massive clusters show 
clear signs of internal chemical enrichment. Further, 
very blue HBs may result from helium enrichment. In 
order to look for signs of helium enrichment in NGC 
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Fig. 11. — Comparison of the color-magnitude diagrams for NGC 2419 and M13 (Gordon et al., in preparation). NGC 2419 data 
come s fr om ACS WFC images, and F435W and F814W magnitudes were transformed to Johnson B and Cousins / using equations 
from [Sirianni ct al. (2005). Only the best measured stars (|SHARP| < 0.5, CHI < 5, and crowding effects totalling less than 0.3 mag) 
have been plotted for NGC 2419. The data for M13 has been windowed to place the blue HB (near the instability strip) level with the 
corres ponding por tion of NGC 2419's HB. Arrows show the positions of possible gaps on the HB, and dotted lines show M5's instability 
strip l ISandQuist fc Bohc 2004), shifted for the higher reddening of NGC 2419 (A(_B - /) = 0.18) . 



2419 using the population ratio R, we need to correct 
for the effect of its blue HB first (in case it has a differ- 
ent cause, such as mass loss) . Helium-enriched stars have 
longer lifetimes than unenriched stars at a given T^g — 
at the blue end of the HB the difference is almost 10% for 
Ay = 0.05. Detection of a helium-enriched subpopula- 
tion of HB stars will still depend on the relative numbers 
of enriched stars, but the correction will provide a more 
accurate idea of the size of such a population. Some au- 
thors have previously att e mpted to correct for the life- 
time effect. ICassisi et al.l (|2003D made corrections to R 
for clusters with blue HB tails by using an HB evolu- 
tionary timescale appropriate to the mean HB mass. For 
clusters with bimodal HB distributions, this can lead to 
significant err or because the lif e time i s not a linear func- 
tion of mass. ICastellani et al.l (j2007l ) attempted to im- 
prove their analysis of uj Cen by applying corrections to 
the HB population using "typical" absolute lifetimes for 
stars in three well-populated parts of the HB. We believe 



our method is a significant improvement in accuracy. 

Although color- Te// relationships remain imperfect, 
it is much more reliable to derive Teff from photome- 
try than it is to derive stellar mass (which is the only 
other plausible choice for a physically-based coordinate 
describing position along the HB). To do this, we used 
DSEP HB models to determine an theoretical "average" 
HB line: specifically the locus of points halfway through 
the HB lives (ending at Yc = 0.05) of stars of different 
masses. The color- Teff transformation we used is shown 
in Fig. [T71 This average line was fit to the stellar pho- 
tometry and used to derive the offset between the theo- 
retical and observed colors. For this cluster, we used the 
FA35W - F8UW color (giving a color shift of 0.13) due 
to its good temperature sensitivity across the color range. 
Note that the evolutionary stage does not have a signif- 
icant effect as long as the color and Teff are determined 
consistently — the other atmospheric variable affecting 
the colors (surface gravity) does not change greatly dur- 
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Fig. 12.— Comparison of the color-magnitude diagrams for NGC 280 8, NGC 2419 a nd M15. NGC 2419 data comes from ACS WFC 
images, and was transforme d to Johnson BV using equations from S irianni et al.l I I2005I ). The same quaUty cuts were applied to the NGC 
2419 photometry as in Fig. [11] with an additional restriction on distance from the cluster center (r < 40"). Photometry for NGC 2808 
(proposal ID 6095; P.I. Djorgovski) and M15 (proposal ID 5324, P.I. Yanny) comes from WFPC2 and was transformed using equations 
within HSTPhot (December 20, 2004 update; sec http://purcell .as.arizona.edu/wfpc2_calib/2004_12_20.htnil). The data for NGC 
2808 and M15 has been windowed to place the blue HB (near the instability strip) level with the corresponding portion of NGC 2419's HB. 



ing the HB phase. 

An additional benefit of using photometry to define 
the temperature scale is that there are features (such 
as the edges of the instability strip) that can be used 
as reference points. For example, the Grundahl it-jump 
has been found in ma ny clusters with blue HB stars 
with Teff > 11500 K (GrundahLeTn] |1999). The u 
jump is difficult to find in NGC 2419 because the HB 
is not heavily populated near its expected position, but 
in M15 it can seen at F336W - F555W = -0.6 (see 
Fig. [TH]). For NGC 2419, the jump appears to be 
at i^336W^ — F555W — —0.65, which corresponds to 
FA39W ~ 21.35 in the WFC photometry, and a small 
color shift can be discerned at this position in the opti- 
cal CMD. This is within the modestly populated blue 
wing of the main peak in the HB distribution. Our 
conversion from F435VF — F81AW independently pro- 
duced Tcff ~ 10950 K, gratifyingly close to 11500 K 
(|Grundahl et~all 119990 . (Tes = 11500 K was predicted 
to be about 0.05 mag fainter.) In addition, the normal 
HB also seems to have a reasonably well- defined termi- 
nation near 29000 K (jBrown et al.il2008f) . As seen in 
Fig. [TI] the end of M13's HB (which appears to corre- 
spond to the end of its hot HB) translates to the middle 



of the EHB clump in NGC 2419 {FA35W « 24). Our 
color-temperature conversion again produces tempera- 
tures that approximately (~ 28500 — 29000 K) match 
the theoretical value. A temperature of 35000 K (corre- 
sponding to the transition in atmospheric helium abun- 
dance spectroscopically observed in NGC 2808 and uj 
Gen) would faU near F814W^ = 24.9 and i^435Vl^ = 24.5. 

As expected from the analysis above, the use of dif- 
ferent sets of models will have a small effect on the 
lifetime-corrected HB star total. For our final analy- 
sis, we fitted a second-order polynomial to the weighting 
factors computed from DSEP models for [Fe/H]= —2.0, 
[a/Fe]= -1-0.2, and Y — 0.244 in order to smooth out 
numerical noise. In addition, we assumed that stars red- 
der than the main peak had evolved from the peak, and 
so were assigned weights from the main peak. We also 
needed to extrapolate the weighting corrections in the 
range 4.4 < logTe// < 4.5, due to a lack of tabulated 
models. This is somewhat risky because a discontinuous 
change in the chemical composition of the atmospheres 
(as there is in NGC 2808 and oj Cen) would cause a dis- 
continuity in the temperature scale derived from the pho- 
tometry. In this case we are likely to be underestimating 
Tcff and the lifetime correction. This would produce a 
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Fig. 13. — The light curve for the bright RR, Lyrae variable V2 
from i?5T photometry from proposal ID 7630. The time coordinate 
is modified Julian date (heliocentric Julian date minus 2450000). 
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Fig. 14. — The effects of blending on the observed properties 
of EHB/BHk stars in the ACS WFC data. A DSEP isochrone 
(jDotter et al. 20070 with age 13 Gyr, [Fe/H]= -2.1, and empirical 
color transformations was fit to the main sequence and lowrer giant 
branch {red line), and the curves show blends of the EHB star with 
a star whose properties were drawn from the fitted isochrone. 

slight bias toward detecting a helium enhancement, but 
our corrections should still remove the lifetime effect to 
first order. With corrections for both incompleteness and 
differences in lifetimes, our original sample of 968 HB 
stars becomes N'j^g = 821 ± 29. (The quoted error bar 
is based on Poisson statistics for the corrected HB num- 
ber.) It is important to note that this would correspond 
to the number of HB stars relative to bright RGB stars 
in our sample if i) the HB stars landed on the ZAHB 
near the blue edge of the instability strip, and ii) the ini- 
tial chemical composition of all stars was the same. For 
this reason, it is appropriate to use this in the helium 



Fig. 15. — Top panel: Detected (solid line) and incompleteness- 
corrected (dotted line) HB star distributions for stars in the ACS 
fields. The numbers refer to incompleteness-corrected numbers of 
stars in each of the magnitude ranges separated by dashed lines. 
Bottom panel: HB star distribution corrected for lifetime differ- 
ences, and including deprojections of evolved stars from the main 
peak {FSliW = 20.1) and probable EHB/MS blends. Inset: Num- 
ber of EHB stars, corrected for incompleteness and blending, and 
including Poisson uncertainties on each bin. 
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Fig. 16. — HB hfetimes as a function of T^f f from DSEP 
Potter et al.l 120071 1 models for (from top to bo ttom) [Fe/H] = 
-1.0,-1.5, -2.0, and -2.5 ([a/Fe] = +0.2), and ISweigar^ ifTOSTl ) 
models Y = 0.30, 0.25, and 0.20 (Z = 0.0001). Le/i panels: lifetimes 
in Myr. Right panels: lifetimes normalized to the stars reaching 
the zero-age HB at logT^ff = 3.85. (For the DSEP model with 
[Fe/H] < —2.0, the most massive models do not reach that tem- 
perature, and were normalized at the lowest effective temperature 
present in the model.) In the top right panel, the bold line shows 
the polynomial fit to the [Fe/H]= —2.0 models that was used in 
making lifetime corrections. 
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Fig. 17. — ACS WFC color-Toff conversion derived from DSEP 
Potter et al.ll2007l ) models for [Fe/H] = -2.0 and [a/Fe] = +0.2. 
The lines correspond to the zero-age HB (dotted line), the HB 
termination (core helium Yc = 0.05; dashed line), and halfway in 
time between the two. No reddening has been applied to the color. 

indicator R (see ^3.5^ . 

If the stars were formed with the same composition, 
then we can use the hfetime-corrected histogram to de- 
rive the production rate of EHB stars relative to other 
HB stars: 33%. However, increased metal and helium 
abundances both increase the HB lifetimes, so any chem- 
ical enrichment among some fraction of the cluster's stars 
would increase N'^^g above the value for a chemically 
uniform population. So the fractional production rate 
of EHB stars in NGC 2419 is at most 33%. Fairly large 
changes in composition are theoretically predicted to pro- 
duce fairly small changes in HB lifetime though: from 
DSEP models, a 0.5 dex chang e in [Fe/H] cha nges the 
hfctime by about 5%, and from lSweigartl ()1987| ) models 
a 0.05 change in Y produces a 7% change among stars 
at the hot end of the HB. Only extreme chemical enrich- 
ment among the EHB stars is likely to mean that fewer 
than 25% of HB stars being produced are being made 
EHB stars. 

We also note that if EHB stars are not produced at a 
constant rate through a cluster's lifetime that an EHB 
population would take some time to equilibrate because 
production rates (relative to other HB stars) would not 
initially equal "destruction" rates (stars evolving off the 
EHB). EHB populations could be expected to linger 
an extra 30 Myr or so after the production mechanism 
(whatever it might be) is shut off. 

3.1.2. An EHB Gap? 

iBrown et al.l (|2001h presented evidence that the faint 
gap in NGC 2808 {B « 20; see Fig. [I2|) corresponds to a 
theoretically-expected gap between EHB stars and blue 
hook stars whic h ignite helium while on the white dwarf 
cooling curve. iMoehler et al.l (l2004l) find that the gap 
in NGC 2808 corresponds to a spectroscopic transition 
from stars with helium-poor atmospheres (brighter than 
the gap) to ones with helium rich atmospheres. They 



find that the brighter stars have temperatures and grav- 
ities consistent with normal HB evolution, and that he- 
lium has probably diffused out of their atmospheres. The 
fainter, helium-rich stars have properties consistent with 
a mixed e nvelope associ ated with a late helium flash. In 
addition, Moehler et al.l found that stars in lo Cen show a 
transition in atmosph eric helium abundanc e at the same 
Teff (~ 35000 K), and lMoehler et all (|1997D found a faint 
HB star in M15 {B = 19.75) with a helium rich atmo- 
sphere, and this star also falls below t h e leve l of the gap 
in NGC 2808. However, iBrown et all (|2001h also noted 
that a gap similar to the one in NGC 2808 does not 
clearly show up in the photometry of some of these same 
clusters, so we have examined this question again here. 

Blue hook candidates have been identifed in several 
clusters through a comparison of CMDs. Photometry in 
the F255W filter appears to give some of the clearest ev- 
idence of BHk stars. Clusters like M13 and M80 have 
HBs that terminate approximately 1.4 magnitudes be- 
low the knee (the blue HB stars with maximum bright- 
ness in F255W ), and thought to ha ve EHB stars but 
no BH k stars (jFerraro et al.|[l998h . iDalessandro et all 
(l2f)f)8h used (F255IF, F255W - U) CMDs to compare 
NGC 6388 and M80, and show that NGC 6388 had a 
group of stars that were fainter than the faintest HB 
stars in M80 and were separated from the rest of the HB 
stars by a gap. Photometry for NGC 2419 in the F250W 
filter from HRC camera (Fig. [3]) places the bright end 
of the EHB about 0.9 mag fainter than the knee and the 
presumptive EHB gap about 1.1 mag below the knee. 
Because the F250W and F255W bandpasses are not iden- 
tical and the photometry for these hot stars are sensitive 
to it, we need a more exact comparison. 

We used the "i? knee" as a referenc e point to shift the 
CMDs of the different clusters, as did 'Daless andro et al.l 
(p008) in their comparison of NGC 2419 with u Cen. The 
knee is essentially based on surface temperature. NGC 
2808's HB is not well populated at the knee, which could 
produce some error, but its red HB allows us to check 
that the alignment with M15 is not far off. The most 
clearly visible gaps in NGC 2419's HB distribution do 
not match up with those in other clusters when the HB 
"knees" are aligned, as in Figs. [11] and [121 The pos- 
sible exception is a marginally visible gap at i? « 24.1 
(J w 24.5) in NGC 2419. Its presence is masked in the 
upper panel of Fig. [15] because of the presence of prob- 
able EHB/MS blends. In the lower panel of that figure 
we have therefore deprojected stars back to the HB fidu- 
cial line using the blending lines shown in Fig. [14] After 
this procedure is executed as part of the incompleteness 
correction process, we see clearer evidence of a gap. The 
statistical significance of such a feature is hard to gauge 
because there is no theoretical expectation for how the 
stars should be distributed, but in the inset we plot cor- 
rected star numbers with Poisson distribution uncertain- 
ties. The bins at F^lAW 24.4 and 24.5 appear to be 
several standard deviations below the general trend de- 
fined by bins on either side. How ever, experiences with 
other clusters (jCatelan et al.ll998l ) lead us to believe that 
there is a significant chance that it could be produced by 
statistical fluctuations even with the relatively large sam- 
ple of stars we have in NGC 2419. However, NGC 2419 
is not the only cluster to show a feature here. 

As shown in Fig. [TT] M13 has produced only a few 
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Fig. 18. — Comparison of the color-magnitude diagrams for NGC 2419 and M15 using the WFPC2 camera and F336W filter. The data 
for M15 has been windowed to place the blue HB (near the instability strip) level with the corresponding portion of NGC 2419's HB. Blue 
crosses are stars that were identified as HB stars via the methods described in ^\ 



possible HB st ars hotter than the clump at the end of its 
blue tail (Ferr aro et al.lfTOQTb . which terminates at the 
level of the weak gap in NGC 2419. In Fig. [H M15 
shows a scattering of stars in the EHB (the one spectro- 
scopically identified BHk candidate falls below the level 
of NGC 2808's gap), uj Gen has a large popula t ion o f 
EHB stars with no gap, although iBrown et al.l ()200lD 
mention that one might have been blurred by the wide 
metallicity dispersi on from star to s tar. S till, a feature 
can be seen in the iGastellani et al.l (|2007( ) HB distribu- 
tion for oj Gen at i? ~ 18.5. The distributions for uj 
G en and NGG 2419 are com pared side-by-side in Fig. 4 
of iDalessandro et all ([2008'), shifted in the same way as 
in Fig. [121 The feature appears there not so much as 
a gap, but as a strong contrast in GMD density: the 
fainter blue hook stars are much more common than the 
slightly brighter EHB stars. Using F-band and shifting 
the GMDs to a ccount for differ ences in distance and red- 
dening, Ro senberg et al.l (|2004f ) found that a gap on the 
extreme blue tail of M54 precisely overlaps the one in 
NGG 2808. The end of the theoretical zero- age HB fit- 
ted to NGG 2419 in Fig. [20] also appears to terminate 
at the approximate position of the apparent gap. All of 



this leads us to postulate that the weak gap at i? ~ 24.1 
in NGG 2419 is real, and marks the boundary between 
the canonical HB and blue hook stars. If true, this fea- 
ture appears to be largely metallicity-independent, and 
its visibility subject to the vagaries of the processes dis- 
tributing stars into the blue tail. 

If the initial chemical compositions of the bluest HB 
stars match those in the main HB peak, then we find 
that 25% of stars reaching the HB are fainter than this 
gap (F8141F > 24.4) and may be blue hook stars. If they 
are more helium-rich, then this fraction would be reduced 
somewhat because we would have undercorrected for star 
lifetimes to get the distribution in the lower panel of Fig. 
1151 The additional correction is very unlikely to change 
the fact that this fraction is considerably larger than in 
any other cluster with the possible exception of uj Gen. 
We will return to this subject in the conclusions. 

3.2. Ultraviolet Bright Stars 

The source of the UV upturn at wavelengths shorter 
than 2300 A in spectra of elliptical galaxies remains 
an unsolved problem in stellar populations. The leading 
candidates for this emission are stellar types that have 
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an appropriate combination of high surface temperature 
and relatively long lifetime. This includes extreme HB 
stars and their evolved forms, the AGB manque stars; 
post-early AGB stars, which evolve to the AGB, but 
leave again before having a thermal pulsation phase; and 
post-AGB stars, which are short-lived but have high lu- 
minosity. The details of the evolution depend largely on 
the m ass of the envelope o utside of the fusion shells. Re- 
cently iBrown et afl (|2008D produced a UV CMD of the 
nearby elliptical galaxy M32, and found that current evo- 
lutionary models predict many more UV-bright post-HB 
stars than are observed. 

NGC 2419 provides an interesting case study because 
of its large mass (which helps to populate relatively short 
evolutionary stages) and because the HB distribution can 
be described to first order with two strong peaks, as dis- 
cussed in WS.1\ Our simple expectation is that the peak 
just to the blue of the instability strip is responsible for 
the majority of the AGB stars observed in the cluster. 
However, an examination of the AGB population may 
reveal whether the AGB stars evolve to a thermal pul- 
sation termination, or whether they leave "early" at a 
lower luminosity level on the AGB. Further, a census of 
stars brighter than the HB and hotter than the AGB can 
help observationally identify how the EHB stars in this 
cluster evolve. 

We have therefore used a combination of optical and 
UV photometry to identify probable hot post-HB stars. 
These are among the brightest objects in the F300W fil- 
ter, as can be seen in Fig. p!9l We first identified the two 
stars cataloged by IZinn et all ([1972) in their search for 
stars that were bright in both the ultraviolet {U, which 
is similar in bandpass to the F300W filter) and visible 
(V). ZNG 2 is more than 30 times brighter than the 
HB in the F300W filter, and ZNG 1 is approximately 1.2 
mag fainter than ZNG 2 and is brighter in the optical. 
Both have optical photometry consistent with being in a 
hot post-AGB phase. We identified five other stars that 
sit closer to the AGB in the CMD. The star PA 4 has 
F300W between ZNG 1 and 2 and is brighter in F555W 
than either ZNG star. Two stars PA 6 and PA 7 did not 
have any UV photometry, but occupy similar positions in 
the optical CMD. Star PA 5 is faintest in the UV because 
it is the coolest star. Star PA 3 is significantly fainter in 
optical bands, and may be an example of a post-early 
AGB star. Star PA 6 was saturated in the WFC images 
(the only field in which it was observed), but the satu- 
rated photometry is accurate enough to conclude that it 
is post AGB. 

There are several reasonable evolutionary scenarios 
that could explain the C MD positions of at least some of 
these stars (see Fig. 4 of lBrown et al.ll2001l for examples 
of evolutionary tracks). As described earlier, an RGB 
star that loses most of its envelope can leave the RGB 
before having a core He flash, evolving rapidly toward 
higher temperatures at nearly constant luminosity. Such 
stars can have a late He flash that lands them on the blue 
end of the HB, or can become a He core WD. Alternately, 
stars that had normal HB and AGB phases can undergo 
H and He shell flashes that lead to large excursions in 
both luminosity and surface temperature. These stars 
would be identified as post-AGB (PAGB) or post-early 
AGB (PEAGB) stars. PAGB and PEAGB stars move 
mostly in Tcff during later flashes, but can get several 



times less luminous than the TRGB for shorter periods 
of time. Post AGB evolution is more likely to explain 
most of these stars, but the possibility exists that one 
of the more luminous stars might be a post RGB star. 
We expect that spectroscopy would the easiest way of 
identifying such stars. 

One of the clearer features associated with the evolu- 
tion of the EHB stars is a group of stars with —0.5 < 
i^435P^ - F8UW < -0.3 in the range 22 < F435M^ < 
23.4 (see Figs. [5] or [T31 for example). We have ver- 
ified that all of the stars to the blue of the inter-gap 
HB population (-0.2 < - F814VF < 0) remain 

to the blue of the HB in all of the UV CMDs. In the 
shortest wavelength filter (F250W), the stars remain a 
nearly vertical sequence, implying luminosity evolution 
from the EHB rather than slight temperature evolution 
from the inter-gap HB (see Fig. [3]). The large size of 
the EHB population compared to the inter-gap popula- 
tion also argues in favor of evolution from the EHB — 
post-HB evolution is quite rapid (on a Kelvin-Hclmholtz 
timescale) in comparison to the HB phase itself, so that 
there are likely to be many fewer post-HB stars. The 
ratio of evolved HB candidates to inter-gap HB stars 
is only about l-to-2, whereas it is about l-to-35 when 
compared to the EHB population. Theoretical ratios of 
post-HB (helium s hell fusion) to HB l ifetimes are typ- 
ically around 0.12 (jCassisi et al.ll2003[) . and it is quite 
possible that not all of the EHB stars produce post HB 
stars in this part of the CMD. DSEP model tracks for 
[Fc/H] = —2.0 and [a/Fe] = -t-0.2 approximately parallel 
the line traced out by these stars in optical bands (see 
Fig. [20]) . Their O.52A/0 track comes close to reproducing 
the star colors, although there are no corresponding HB 
stars on the slowest portions of the track (between the 
zero- age HB and the portion shooting toward higher lu- 
minosity after core helium exhaustion) . The gap of about 
0.5 mag in the F^ZhW filter is expected to correspond 
to the rapid evolution between core helium exhaustion 
and the establishment of a stable He fusion shell. This 
may be useful as an envelope helium abundance indi- 
cator specifically for the EHB stars because the helimn 
abundance affects the envelope mass for a given effective 
temperature, and the envelope mass affects the star's lu- 
mino sity when the He shell is established (jBrown et al.l 
[20081) . For more massive envelopes, the H fusion shell 
can be restarted at the same time, and this contributes 
significantly to the luminosity. At the moment, publicly- 
available models for EHB stars (evolved past central He 
exhaustion) with helium enrichment do not exist, so we 
are unable to calibrate this. We encourage such model- 
ing because it may help answer questions related to the 
origins of multiple populations in clusters and the possi- 
bility that it is caused by cluster self-enrichment. 

A second group of stars with F555W^ w 19.5 also ap- 
pears to be composed of evolved HB stars. These appear 
to correspond to "supra-H B" stars identified in cl usters 
like M 13 and discussed in IStrom et al.l (|1970f ) and 'Zimi 
(Il974l ). The possibility that most of them are optical 
blends can be ruled out based on their brightness in the 
ultraviolet. (One star that falls in this range was identi- 
fied as a blend in the HST photometry, however.) Among 
these stars is the bright RR Lyrae star V2, which was 
discussed in ^3.1] We consulted theor etical HB models 
(jDotter et al.ll2007t ICassisi et aIll2004D with very metal- 
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Fig. 19. — CMDs of ultraviolet-bright post-HB stars. The triangular symbols point in the expected direction of the evolution: upward 
for AGB manque stars, rightward for supra-HB stars evolving toward the AGB, and leftward for post-AGB and post-EAGB stars. 



poor compositions and no helium enrichment, and veri- 
fied that stars originating from the primary peak in the 
HB distribution are expected to produce AGB stars ex- 
clusively — the evolution of the stars keeps them within 
about 0.1 mag of the HB until central He exhaustion 
when they move to the AGB clump in about a Kelvin- 
Helmholtz timescale. Models agree that as stellar mass 
decreases the post-HB evolution toward the AGB occurs 
brighter and brighter in V until the tracks never reach 
all the way to the canonical AGB. In a small range of 
masses, stars can eventually evolve most of the way to 
the AGB, but spend a significant amount of time (more 
than a Kelvin-Helmholtz timescale) brighter than the HB 
at intermediate temperatures before reaching the AGB. 
Such a phase would correspond to the AGB clump if the 
stars had reached the canonical AGB. It is also related 
to the "bump" on the first ascent giant branch — during 
this time the He fusion shell comes into equilibrium af- 
ter consuming the He abundance gradient left by the HB 
phase and adjusting to the new He fuel concentration. 

A big question is where on the HB the supra-HB stars 
originate from. According to canonical models, the only 
evolution tracks that pass through this region come from 
inter-gap HB stars. However, the relative numbers of 
inter-gap stars and supra-HB stars appears to be incon- 
sistent with this. Fig. shows a comparison between 
our ACS WFC photometry and evolutionary tracks and 



synthetic HB populations from DSEP models. [We used 
DSEP products using their "empirical" color transfor- 
mations as these appear to do a better job than the 
"synthetic" transformations in repr oducing ACS data 
in and around the su bgiant branch (iDotter et al.l [20071 : 
ISaraiedini et al.|[2007[ ) .] The numbers of stars in the syn- 
thetic model representing the inter-gap population does 
not generate enough supra-HB stars (especially on the 
hot side), although it may explain the reddest stars in 
that group. Our best explanation for the number of 
evolved stars that are seen is that both the supra-HB and 
AGB manque stars have evolved from the EHB peak, 
with the supra-HB evolving from the cooler (brighter in 
optical filters) part of the peak. 

Increased envelope helium abundance tends to give 
stars a more pronounced blue loop during their early HB 
evolution, and they tend to retain higher surface temper- 
atu res during their subsequent evolution. (For examples, 
see lSweigartlll987l ) Helium enrichment can also lead to 
a conver gence of post -HB evolutionary tracks (see Fig. 
9 of Buss o et al.ll2007j for an example) that could be re- 
sponsible for supra-HB stars. For NGC 2419, an increase 
in envelope helium among the EHB stars might bring 
the models and observations into better agreement, but 
the increase may not need to be large since the theoret- 
ical tracks come close to matching the properties of the 
evolved HB stars. Once again though, theoretical models 
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Fig. 20. — Comparison of ACS WFC photometry with DSEP HB models. In all panels theretical values are shifted by 19.97 in F555W and 
0.2 in F439W-F814 to account for distance and reddening. Left panels: Comparison with evolution tracks having [Fe/H]= —2.0 and [o/Fe] 
= +0.2. The ZAHB is shown in red. The tracks are for 0.51, 0.52, 0.54, and 0.56 A/q {bottom panel), 0.58 and 0.6 Mq (middle panel), and 
0.65 and 0.7 Mq {top panel). Right panels: Examples of synthetic horizontal branch populations having [Fe/H] = —2.1, [a/Fe] = +0.2, 
and a mass dispersion of 0.01 Mq, The number of synthetic HB stars roughly corresponds to the number in the incompleteness-corrected 
sample (see Fig. I15|l . Also included [solid line) is an isochrone for the same composition and an age of 13 Gyr. Crosses show our mode fits 
to MS data, and individual stars that represent the mean HB line. 
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for EHB stars with varying helium abundance are need 
to settle the question. 

3.3. The Upper Red Giant Branch 

Because NGC 2419 has a significant population of 
EHB stars, it is natural to ask why these stars have 
been produced in such large numbers. In recent years, 
there has been growing evidence that subpopulations of 
chemically-enriched stars have been responsible for mul- 
timodal distributions of HB stars and for the existence of 
stars at the high temperature end of the HB. However, 
there has been no evidence of chemical composition vari- 
ations within NGC 2419 from the CMDs (no large disper- 
sion along the evolutionary sequences as seen in uj Cen, 
and no sign of multiple sequences), although it should 
be said, there has also been little in the way of spectro- 
scopic abundance studies. As a result, we examined the 
possibility that extreme mass loss could be responsible 
for the bluest stars in the cluster. Helium enrichment 
within a stellar population does encourage the formation 
of blue HB stars by allowing lower mass stars to evolve 
off the main sequence and up the RGB. However, helium 
enrichment does not change the theoretical requirement 
that nearly all of the mass outside the helium core must 
be removed to put it on the EHB. If bright RGB stars 
deplete the mass stored in their envelopes before reach- 
ing the TRGB, they can avoid a traditional helium flash 
event, and instead may have a "hot flash" helium igni- 
ti on that ultimately lands them on the EHB (see Fig. 4 
of lBrown et al.l[200lL for example) . Stars that can some- 
how lose their envelopes before reaching within about 0.4 
magnitudes of the TRGB will not subsequently ignite 
helium, and are expected to produce He white dwarfs 
(jD'Cruz et al.l[l99fiHBnwn et al.l[2?)nl . 

Infrared excesses associated with giants in globular 
clusters provide evidence of dust fo rmation, poss i bly as - 
sociated with strong mass loss. lOriglia et al.l (|2007l ) 
found excesses among a fraction of stars in the metal- 
rich clusters 47 Tuc at luminosity levels down to the HB, 
and becoming stronger with luminosity. The presence 
of excesses among only a fr action of th e stars implies 
episodic mass loss, although iBover et al.l (|2008[ ) present 
evidence that much of this may have been due to blending 
in the infrared photometry. iBover et al.l find no evidence 
of mass loss below the TRGB in ui Cen, which is a cluster 
that is closer in composition to NGC 2419 that also has 
a large population of EHB stars. So at this point there is 
not strong evidence of strong mass loss below the TRGB. 

We used the cumulative luminosity function in our 
analysis, counting stars starting at the observed TRGB. 
Comparison with models requires a magnitude shift and 
a vertical normalization. The horizontal magnitude shift 
was taken from the magnitude difference between the 
brightest observed red giant of the cluster and the TRGB 
in the theoretical model. By doing this, we virtually elim- 
inate all depende nces on age and che mical composition 
for [Fe/H] < -1 ( Sandquist fc Martel [20071 ). It is partly 
for these reasons that the luminosity of the TRGB in / 
has been used as a distance indicator for large, old metal- 
poor stellar populations. While there are giant stars that 
are brighter than the TRGB we identify, there is a large 
gap in FSlAW that leads us to identify them as AGB 
stars, as shown in Fig. [^TJ The CMD also appears to be 
more uniformly populated fainter than the third bright- 
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Fig. 21.— (F8UW, F555W - FSliW) CMD for NGC 2419 
stars on the bright RGB and AGB. The photometry is averages of 
measurements on the HST flight system after individual datasets 
were corrected for zeropoint differences relative to the proposal 
7630 data. 

est giant star having F814W ^ 15.78 that we identify 
with the TRGB. 

For the vertical normalization we forced the models 
to have the same number of observed stars just brighter 
than the RGB bump at F814W^ = 18.4. Shown in Fig. 
[22l are the observed data matched with models having 
[Fc/H] of -1.84 and -2.01. Although the [Fe/H] = 
— 1.84 model does the best job in matching the position 
of the RGB bump relative to the TRGB (see the inset in 
Fig. [22), in both cases the models predict a larger num- 
ber of red giants than the observed data has between 
F8UW 16.4 and f814 W^ ~ 17.7. As discussed by 
ISandauist fc MarteH (|2007f ). the shape of the cumulative 
luminosity function is a robust theoretical prediction. 

The models were compared using the Kolmogorov- 
Smirnov (K-S) test. This evaluates the statistical prob- 
ability that the null hypothesis (the two distributions 
are drawn from the same distribution function) can be 
disproven. The probability is based on the D statistic, 
the maximum absolute difference between the cumulative 
distributions of the two sets. Although it has the advan- 
tage of making no assumptions about the precise form 
of the distributions, the probability that is calculated 
depends on the size of the sample. Because the choice 
of the cutoff magnitude for the comparison changes the 
sample size, we ran the K-S test with different cutoffs (es- 
sentially different normalization points) . Table [3| shows 
the results of the K-S tests, including the number of stars 
brighter than each magnitude cutoff. 

NGC 2419 only shows a slight deficiency of giants com- 
pared to the theoretical models. The minimum K-S test 
probability is 16% for a sample with F8UW « 17.9. The 
maximum difference between the observations and mod- 
els (0.09 in logio N) seen at i^814W^ = 16.6 in Fig. [22| 
is comparable to the fraction of EHB stars relative to all 
HB stars; a 33% decrease in the giant population would 
only correspond to a decrease in logj^Q of about 0.12. 
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Fig. 22. — The cumulative luminosity function for bright 
RGB stars for NGC 2419, along with Victoria-Regina models 
(|VandenBcrg ct al. 2006;) for age 12 Gyr and [Fe/H] = -1.84 (thin 
solid line) and —2.01 (thin dotted line) The inset shows a zoom 
centered on the position of the RGB bump. 



TABLE 3 

Cumulative Luminosity Function Comparisons 
FOR, NGC 2419 



AFSUW 


N 


D 


P 


D 


P 






[Fe/H] 


= -1.84 


[Fe/H] 


= -2.01 


1.0 


41 


0.102 


0.769 


0.097 


0.817 


1.2 


54 


0.067 


0.962 


0.068 


0.958 


1.4 


72 


0.101 


0.434 


0.100 


0.441 


1.6 


85 


0.060 


0.911 


0.058 


0.930 


1.8 


102 


0.048 


0.970 


0.046 


0.981 


2.0 


131 


0.073 


0.480 


0.072 


0.492 


2.2 


154 


0.055 


0.721 


0.054 


0.744 


2.4 


181 


0.043 


0.873 


0.045 


0.852 



Note. — Col. 1: magnitudes below TRGB for cutoff 
of RGB sample. Col. 2: number of RGB stars in sam- 
ple. Cols. 3, 5: absolute deviation between observed 
and theoretical cumulative luminosity functions. Cols. 
4, 6: Probability from K-S test. 



As a whole though, the resuhs are marginally consistent 
with a population of giant stars having a distribution the 
same as the theoretical luminosity function. The results 
appear to validate the algorithms for computing plasma 
neutrino losses near the TRGB — Sandquist & Martel 
([2007) no ted that ne wer algorithms (e.g. iHaft et al.1 
[1994: ,Itoh et al.lll996D produce LF shapes that are con- 
sistent with observations. For typical RGB core condi- 
tions the rapid increase in neutrino energy loss is strongly 
dependent on the large temperatures and densities, and 
neutrino emission rapidly increases as the TRGB is ap- 
proached. The shape of the LF would be changed if 
there were any energy loss mechanism having a signifi- 
cantly different dependence on temperature and density. 
The good agreement of the present NGC 2419 data with 
the models strengthens the case that the plasma neu- 
trino losses are now being accurately modeled and non- 
standard energy loss mechanisms are not operating in 



the core in a significant way. 

The position of the modest deviation in NGC 2419's lu- 
minosity function also shows up at a different luminosity 
than t he one noted in NGC 2808 bv lSandauist fc IVIartell 
()2007D . As noted earlier, stars within brightest 0.4 mag 
of the luminosity function are expected to ultimately pro- 
duce HB stars: RGB stars having a normal helium flash 
at the TRGB will become normal HB stars, while RGB 
stars leaving the RGB a little early can have a late flash 
and become EHB stars. In NGC 2808, the brightest 
1 mag in F814W seemed to be fairly uniformly below 
predictions by about 20%. In NGC 2419 the brightest 
0.6 mag seems to be in good agreement with the mod- 
els, while there are noticeable deviations in the next 1.4 
mag. Even taking into account the longer lifetimes of 
blue HB stars having enhanced metal or helium abun- 
dances, the HB distribution implies that more than 25% 
of all HB stars being produced in NGC 2419 are becom- 
ing EHB stars (see gSTj). In NGC 2808, the fraction of 
HB stars that are on the EHB is signiflcant, but quite a 
bit smaller (see Fig. [T^l for example). Based on a simple 
interpretation of the RGB luminosity function, we would 
have expected NGC 2808 to have a large population of 
EHB stars, and NGC 2419 to have a relatively unpopu- 
lated EHB but a potentially substantial population of He 
white dwarfs. On the face of it, these results are some- 
what contradictory. What could be happening in these 
clusters? 

A first possibility is that there is unaccounted-for con- 
tamination of the RGB sample by AGB stars that masks 
a larger deficit of RGB stars in NGC 2419. EHB stars 
are not expected to evolve into stars that would be con- 
fused with RGB stars because they remain at high T^b 
in their post-HB evolution, but the main group of HB 
stars does produce AGB stars falling near the RGB. The 
bright RGB is particularly difficult to decontaminate be- 
cause the AGB approaches the RGB most closely near 
the TRGB, and this is also the part of the RGB where 
contamination would have the largest effect on our ability 
to detect signs of RGB stars leaving early. More accurate 
photometry or spectroscopy would help to separate the 
AGB and RGB samples more definitively. However, our 
AGB star sample is quite large in comparison to the total 
HB population (i?2 = Nagb/N'jj^ = 0.146), and espe- 
cially in comparison to the HB population in the main 
peak. So we believe we have identified the great majority 
of the AGB stars. 

If the luminosity function shape is correct, is it possible 
that production of EHB stars does not require extreme 
mass loss? The primary suspect for accomplishing this 
is through a cluster's self-enrichment of helium. There is 
quite a bit of evidence supporting the picture of chem- 
ical self-enrichment within clusters, such as the sloping 
HB a nd unusually long period s of RR Lyraes in NGC 
6441 (jCaloi fc D'Antonall2007[ ). There is also some evi- 
dence that chemical self-enrichment can reduce the need 
for dispersion in red giant mass loss to produce observed 
HB color distributions. For example, in M3 (which ad- 
mitted ly does not have EHB stars), ICaloi fc D'Antonal 
(|2008f ) found that the combination of substantial color 
spread on the HB and a sharply peaked period distribu- 
tion for RR Lyrae stars is most consistent with dispersion 
in helium abundance and virtually no dispersion in mass 
loss (< O.OO3M0). However, unless there is an extreme 
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form of helium enrichment in the cluster, it would still 
be necessary for stars to lose nearly all of their envelopes 
to reach the EHB. An initial helium abundance Y ^ 0.4 
would be required to reduc e the mass of stars at th e 
TRGB to - O.5M0 (e.g. Ncws ham fc Terndrupl [2001 . 
An enrichment of this size should produce a clear second 
sequence in the CMD (remember that such a sequence 
would also need to have t o be about 25% of the entire 
population). In NGC 2808. lPiotto et all ()2007l ) identified 
a triple main sequence that was consistent with popula- 
tions having different helium abundances ranging from 
primordial {Y « 0.248) up to 0.35 - 0.4. uj Cen shows 
a bluer main sequence (|Bedin et al.ll2003 ) that is more 
metal-rich than the majority of stars, and there fore is 
probably also helium enriched (jPiotto et al.ll2005[ ). 

If similar helium variations were responsible for the two 
HE populations in NGC 2419, a bimodal main sequence 
should be present in NGC 2419. In Fig. [23l we show a 
zoom on the main sequence in a high-quality subset of 
the ACS WFC photometry. A fiducial line was derived 
from stars in 0.1 mag bins, and the main sequence was 
rectified by subtracting the color of the fiducial line inter- 
polated to the magnitude of each star. Histograms of the 
MS star samples in 0.25 mag bins are shown on the right 
side of the figure. The distributions resemble Gaussian 
distributions, and no feature appears consistently from 
magnitude bin to bin. A large chemically-distinct popu- 
lation can be ruled out with the present data unless there 
is a dysfunctional combination of helium and metal en- 
richment that cancels any color shift. 

Other factors affect the interpretation of present-day 
samples of RGB and HE stars. For example, unless 
there is a physical process creating a semi-permanent 
feature in the luminosity function, an observed deficit 
of RGE stars (say associated with a statistical fluctua- 
tion) would propagate up the RGE, and eventually re- 
sult in a reduction in the number of HE stars, although 
it should have no bearing on the issue of the production 
of EHB stars. The current sample of HB stars comes 
from stars that were on the RGE within about the last 
130 Myr. According to models, a star moves through 
the uppermost 4.5 / mag of the giant branch (almost 4 
mag in V) in 100 Myr, which is a typical lifetime for HB 
stars near the instability strip. For a population with 
a uniform composition and no mass loss, the stars in 
this range only vary in initial mass by about 0.002Mq 
— a tiny range when translated back to the main se- 
quence. In other words, the rapidity of RGE evolution 
should stretch out statistical fluctuations on the MS to 
the point of undetectability. The RGB cumulative lumi- 
nosity function is also insensitive to changes in the total 
number of RGE stars in this range due to the way it 
is normalized. What Fig. [22] examines is variations in 
the numbers of stars at different luminosity levels within 
the sample. It is interesting that the two massive clus- 
ters with EHB stars that have been studied in this way 
have both shown small deficits toward the bright end of 
the RGE (and not excesses over theoretical predictions). 
However, the observed deficits appear at different lumi- 
nosity levels (implying that different kinds of stars are 
being produced) and have m arginal significance. In ad- 
dition, only one cluster (M5: ISandauist fc Marte]| |2007f l 
without EHB stars has been discussed in a similar man- 
ner, although it was consistent with theoretical expec- 



tations. We are forced to conclude that the different 
methods of examining RGB mass loss (detection of cir- 
cumstellar dust and the luminosity function of the upper 
RGE) do not provide convincing evidence of strong mass 
loss well below the TRGB as yet. 

3.4. Radial Distributions 

ICastellani et al.l ()2007f ) found trends in the radial dis- 
tributions of stars in the massive cluster lu Cen, including 
a radial variation in the percentage of stars occupying the 
blue end of the HE. Structurally NGC 2419 and uj Cen 
are similar in co ncentration (c = 1.67 for NGC 2419 , and 
1.34 for UJ Cen; [McLaug hlin fc van der Marell [20051 ). al- 
though NGC 2419 has unusually large core and half-light 
radii (r^ = 8.7 pc and rh = 19.9 pc) when compared to 
other clusters of similar luminosity (4.1 pc and 7.2 pc 
for LU Cen, for example). However, when compared with 
other clusters of similar galactocentric distance, its core 
and half-light radii seem typical, but its total luminosity 
does not. We have a limited radial extent to examine 
here, but we have compared the cumulative radial dis- 
tributions for the core sample in Fig. [24l A two-sample 
Kolmogorov-Smirnov test returns a probability of 55% 
that the RGB and bright HE stars are drawn from the 
same distribution, and of 47% for the RGE and AGB 
samples. So the RGB, AGB, and bright HE distributions 
are in reasonable agreement, with some small deviations 
of low significance. 

3.5. The R Population Ratio 

Because the population ratio R — Nhb/J^rgb is a 
sensitive indicator of helium abundance for clusters with 
large samples of evolved stars, we will attempt an ad- 
ditional test of the hypothesis that unusual CMD mor- 
phologies in the most massive clusters can be attributed 
to chemical self-enrichment. 

The main difficulty in using the R ratio is the calibra- 
tion of the values. Recent models ([Salaris et al.l 1200^ 
predict R should have a value near 1.40 over a wide range 
of metallicities, but that the precise value depends on 
the HB morphology. To make a reliable comparison of 
He abundances in a relative sense, we must correct for 
variations in lifetimes along the HB. This correction is 
important for a cluster with a blue HB morphology like 
NGC 2419, and we discussed our method in §3.11 

In order to compare NGC 2419 to other clusters, 
we follow a method s imilar to that described by 
iRecio-Blanco et all (|2005f ) to select the faint limit for the 
RGE star sample. In their method, a template cluster 
having an intermediate HE morphology is fit to the clus- 
ter data and used to identify a reference HB level. Be- 
cause NGC_2419is_a^ery metal-poor cluster, we used 
M68 ([Piotto et al.ll2002[ ) as the template. The M68 pho- 
tometry was shifted in magnitude and color to achieve 
the best possible fit. (As discussed in §3.11 and below, 
the red HE stars in NGC 2419 seem to exclusively be 
evolved HE stars, so we looked for the best fits with 
stars in the main HE peak and bluer.) The comparison 
between the proposal 7628 dataset (the o nly one that 
used t he F439W and F555W filters) and the lPiotto et"aLl 
(|2002f l dataset for M68 is shown in Fig. [25] with the 
M68 photometry was shifted by 4.75 mag in F555W 
and 0.05 mag in FA39W ~ F555M^. The error in the 
magnitude comparison appears to be about the same 
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F435W-F814W 



Fig. 23.— Left panel: The NGC 2419 main sequence in the WFC 
Right panels: Histograms of stars from the left panel after the color 
into 0.25 mag subsets. 

as the 0.07 mag quoted by iRecio-Blanco et aD ()2005l ) 
in their work. Based on the reference ZAHB value fo r 
M68 {F^d^W^"^"^ = 15.73: IRecio-Blanco et al.l[2Q0l . 
this gives F555W^^"^ = 20.48 ± 0.08 for NGC 2419. 
Taking into account the zeropoint difference between the 
datasets for proposals 7628 and 7630 (—0.017 mag in 
F555W), we have FbbbW^"^"^ ^ 20.46 ± 0.08 on the 
reference system. The ZAHB level in F555W is used in 
the defini tion of the faint end of th e RGB sampl e in th e 
studies of IZoccah et~aLl (|2000D and lSalaris et all (|2004D . 
and we have used the same definition here for the most 
direct comparison to other clusters. 

With the zeropoint, we find Nrgb = 529^32 for stars 
in the same fields as the FEB stars (i.e. the ACS WFC 
fields of view). The quoted error bars on Nugb come 
from the uncertainty in the ZAHB level, which dominates 
Poisson statistics as a source of error. Although NGC 
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ACS dataset for stars with quahty and radial position (r < 90") cuts, 
the fiducial line has been subtracted. The total sample is separated 

2419 has a fairly large number of RR Lyrae stars that 
could be used to establish the ZAHB level, indications are 
that all of the RR Lyrae stars have significantly evolved 
away from the blue HB (see i |3T|) . iFLipepi et all ()2007l ) 
find (Vrr) = 20.31 ±0.01 from 67 stars, but also find that 
the period-amplitude distribution is most consistent with 
the "well-evolved" RR Lyrae varia bles in M3. For the RR 
Lyrae stars in M68. 1 Walked (fl99l finds (Vrr) = 15.64 ± 
0.01. The M68 RR Lyrae stars have positions in the 
period-amplitude diagram that are consistent with the 
"regular", relatively unevolved stars in M3. The implied 
magnitude shift between NGC 2419 and M68 from (Vrr) 
measurements is 4.67 mag. This is smaller than the value 
derived from CMD comparison, and does not produce an 
adequate fit. This again is consistent with the idea that 
the NGC 2419 RR Lyrae stars are evolved and do not 
provide a good estimate of the ZAHB position, so we 
therefore rely on the CMD fit value. 
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Fig. 24. — The cumulative radial distributions for evolved stars in the combined HST fields. The red line shows bright RGB stars 
(F814VK < 20), green shows AGB stars, and blue shows HB stars excluding the EHB {FS14W < 23.75) 

2419 stars is higher than primordial by about 0.016, but 
the errors indicate this is just above a la difference. This 
is in spite of at least two effects that probably produce 
systematic biases toward larger helium abundances (see 
H3.1^ . (Without the correction for star lifetimes, the HB 
population would have been Nhb — 954, giving R — 1.81 
and implying an average helium abundance enhancement 
of about 0.041.) Still, if the EHB population (about 31% 
of the corrected total) is the only part that is helium en- 
riched, it could be enhanced by about 0.05. Based on 
the evolved stars, there is only weak evidence of helium 
enrichment. 

As a final comparison, we mention the populatio n 
ratios for lu Cen computed bv ICastellani et alT (|2007f) . 
Though they did not compute an R value that could be 
directly compared with ours, their comparisons of HB 
to RGB or MS stars implied that the HB star sample 
was 30-40% larger than predicted by canonical evolu- 
tion models. Their result was actually larger than would 
be expected given that to Cen appears to host helium- 
enriched populations. Regardless of the reason for the 
discrepancy, it further underlines the weak evidence of 
helium enrichment in NGC 2419. 

4. CONCLUSIONS 

We have examined archival HST imagery in the opti- 
cal and ultraviolet to study evolved stars in the core of 
the massive globular cluster NGC 2419. An important 
aspect of these populations is the clear separation of the 
horizontal branch population into two groups: the pri- 
mary peak just to the blue of the instability strip and a 
secondary peak of comparable size at the extreme end of 
the HB. More than 38% of the HB stars in this cluster 
are in the EHB peak, and if their initial compositions 
were the same, this means that more than 33% of its 
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F439W - F555W 

Fig. 25.— WFPC2 CM Ds for NGC 2419 (dots) and M68 {open 
circles: [Piotto et al.ll2003) . The M68 dataset has been shifted to 
match the blue HB of NGC 2419. 

When combined with N'^b = 821 ± 29 (the HB sam- 
ple when corrected for incompleteness and for lifetimes 
of blue HB stars relative to stars at the edge of the insta- 
bility strip), we find R — 1.55^.^'^^. The error quoted for 
N'fjg only includes Poisson errors, but the uncertainty 
in B deriving from the Z AHB magnitude level dom- 
inates the uncertainty in R. For [Fe/H] = —2.1 and 
Y = 0.245. ISalaris et all (|200l find R = 1.40 from theo- 
retical models, and this value changes by at most 0.01 at 
[Fe/Hl= -1.8. This implies that the average Y for NGC 
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red giants have been diverted into this population (more 
than 25% into blue hook candidates). While our study 
only examined the core population (and in uj Cen, for 
example, there is some evidence for a decrease in the 
abundance of EHB stars with distance from the clus- 
ter center), this is a feat that is even remotely matched 
by only some of the very most massive clusters in the 
Galaxy. For example, lo Cen has approximately 2 7% of 
its HB stars on the EHB ()Dalessandro et al.l [20081) — a 
slightly smaller proportion than NGC 2419. M54 also 
has a modest popu lation of EHB stars (|Rosenberg et al.l 
l2004HSiegel et al.l l20G7). but it is a significantly smaller 
fraction of the to tal HB population (we estimate < 15% 
from the CMD in lRosenberg et al.|[200l . NGC 2808 also 
has a modest EHB population that is 12% of its total HB 
sample ( Castellani et al.l 120061 where we have included 
their EBT3 stars and the HBp stars that are likely to be 
composed of blends of EBT3 a nd MS stars). In the ver y 
metal-rich cluster NGC 6388 (jPalessandro et al.l 120081 ) . 
only about 3.5% of the HB stars are on the EHB (2% 
as blue hook stars). A few other massive clusters (e.g. 
NGC 6441, M15) have a handful of EHB stars that com- 
pose a few percent of their populations at most, and an- 
other massive cluster (47 Tuc) has no known EHB stars. 
w Cen, M54, NGC 6388, NGC 6441, and 47 Tuc are 
probably th e only clusters that are more massive than 
NGC 2419 ()McLaughlin fc van der Marel [2005f) . while 
NGC 2808 appears to have a similar mass and M15 is 
approximately 50% less massive. So while large clus- 
ter mass certainly makes the appearance of EHB stars 
more likely, another factor is influencing their relative 
abundance. This factor is probably not past cluster mass 
— NGC 2419 is unlikely to have lost as much mass in 
its history as clusters like uj Cen and M54 because its 
orbit keeps it far out in the Galactic halo, minimizing 
the effects of Galactic tides and disk shocking. There is 
also no evidence of populations of greatly different age in 
NGC 2419 based on a superficial examination — NGC 
1851 shows two distinct subgiant branches that could re- 
sult f rom an age difference of about 1 Gyr (jMilone et al.l 
[200l) . 

For the clusters above, there is a hint that the frac- 
tion of EHB stars anticorrelates with metallicity. M15 is 
about 0.1 dex more metal poor than NGC 2419, but also 
by far the least massive of these clusters, lo Cen has a 
large metalhcity spread (-2.2 < [Fe/H] < -0-7), bu t the 

most common [Fe/H] 1.75 ( Johnson et al.ll2008 '). Of 

the clusters with moderate EHB populations, M54 has 

[Fe/H] 1.6 (iBrown et al. 199S) and NGC 2808 has 

[Fe/H] 1.15 (|Carrettaii2006) . The remainder of the 

massive clusters (NGC 6388, NGC 6441, and 47 Tuc) are 
very metal rich ([Fe/H] > —0.7). This crude correlation 
may be related to well-known relation between [Fe/H] 
and HB morphology — more metal-poor clusters tend to 
have bluer HB stars. 



Even if we compare NGC 2419 to its nearest rival in 
producing EHB stars (w Cen), we find a lack of a clear 
lead. bJ Cen has an HB morphology similar to NGC 
2419, though Lo Cen has diverted a greater proportion of 
stars into the blue HB tail between the two peaks, which 
would imply that multiple populations should overlap 
to a greater degree in uj Cen. w Cen shows very clear 
signs of multiple populations on the giant branch and 
main sequence (Be din et al.ll200'i ) , including a blue main 
sequence that is more metal-rich than the majority of 
stars, and therefo re is probably also helium enriched 
(jPiotto et al.ll2005f l. The number ratio of chemically en- 
riched group to the total MS star population is similar 
to the r atio of EHB s tars to the total HB population, 
and thus iPiotto et al.l suggested that the EHB stars are 
the descendants of the enriched population. In NGC 
2419, there is no evidence to date of multiple popula- 
tions among red giants or main sequence stars. Although 
there is a little evidence consistent with slight helium en- 
richment in the cluster (based on post-HB evolution and 
the R population ratio), there is no evidence that it is 
as strong as in uj Cen. It is potentially possible to hide 
multiple populations on the main sequence by increas- 
ing metallicity and helium abundance in a way to cancel 
out their opposing color shifts. If this is happening in 
NGC 2419, it would imply a smaller helium enrichment 
(Ay ~ 0.03) than is implied for uj Cen (Ay = 0.14). 
Such a small helium enrichment would be consistent with 
the indicators presented in this paper, but would be un- 
able to explain the two HB star populations. Spectro- 
scopic observations of cluster giants would still be needed 
to determine whether a metallicity spread is truly present 
in the cluster. 

In one sense NGC 2419 is like oj Cen: neither ap- 
pears to have undergon e strong dynamical relaxation 
(Daless andro e t al. 200.^ iFerraro et alll20060 . There is 
littl e or no evidence of a correlation with concentration c 
(McLaughlin &: van der Marell l2005l ) though. The three 
massive clusters with the largest EHB fractions (NGC 
2419, Lo Cen, and M54) are also largest in Rh, although 
47 Tuc has similar structural properties to M54 and has 
no EHB stars. 

At this point, we must confess that while NGC 2419 
seems to have less complicated stellar populations overall 
than clusters like lo Cen and NGC 2808, its relative sim- 
plicity helps to stymie attempts to draw out a coherent 
picture of the source of multiple populations and extreme 
horizontal branch stars. 
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TABLE 4 
Evolved Stars in NGC 2419 



ID 


5R.A. 


<5DEC 


F300W 


F555W 


F814W 


F250W 


F336W 


F380W 


F435W 


F439W 


F450W 


F675W 


Notes 




(arcsec) 


(arcsec) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


















Post-AGB candidates 










PA 1 


7.800 


16.468 


16.876 


17.461 


17.371 








17.403 








ZNG2 


PA 2 


30.141 


-28.501 


18.026 


16.632 


15.707 




17.390 


17.617 




17.397 


17.137 




ZNGl 


PA 3 


26.780 


-14.389 


19.339 


17.580 


16.632 


















PA 4 


-49.254 


-100.918 


17.556 


15.789 


14.797 




16.576 


16.851 


16.426 


16.767 








PA 5 


-8.697 


-53.664 


19.832 


16.662 


15.200 




18.842 


18.536 


17.791 


17.921 


17.593 






PA 6 


76.915 


-185.451 














17.594 










PA 7 


-92.211 


-73.163 




16.563 


15.442 








17.828 
























supra-HB stars 










SH 1 


32.980 


-41.293 


18.132 


19.545 


19.653 








19.414 










SH 2 


-33.302 


13.065 


18.755 


19.579 


19.600 








19.574 










SH 3 


10.784 


13.315 


19.941 


19.209 


18.525 








19.661 










SH 4 


-1.433 


-21.805 


20.141 


19.582 


19.374 


20.164 


19.961 


19.862 


19.716 


19.749 


19.763 


19.535 




SH 5 


5.825 


-5.158 




19.411 


19.013 








19.727 










SH 6 


5.765 


-5.179 


19.849 


19.062 


18.259 




19.450 


19.439 


20.146 


19.310 


19.217 


18.432 




SH 7 


-7.038 


-23.101 


20.605 


19.452 


18.748 


21.439 


20.046 


20.272 


19.803 


20.190 


20.077 


19.206 


V2 
















AGB manque candidates 










AM 1 


-63.325 


-13.723 




22.200 


22.490 








22.058 










AM 2 


11.255 


-1.984 


20.762 


22.363 


22.590 








22.217 










AM 3 


-20.198 


-14.353 


20.645 


22.326 


22.501 




20.709 


21.975 


22.217 


22.310 


22.336 


22.466 




AM 4 


25.687 


-57.545 


20.390 


22.599 


22.796 








22.469 










AM 5 


52.320 


-55.610 


20.649 


22.819 


23.059 








22.619 










AM 6 


43.226 


-65.388 


20.914 


23.055 


23.327 








22.884 










AM 7 


5.899 


-27.397 




23.163 


23.403 








22.981 










AM 8 


1.987 


-3.072 


20.954 


23.010 


22.860 


20.551 


21.101 


22.455 


23.433 


22.787 


22.601 


22.445 


blends with BBS in WFPC2 


AM 9 


-40.218 


-0.963 


21.156 


23.312 


23.547 




21.271 


22.809 


23.126 


23.184 


23.130 


23.400 




AM 10 


-32.261 


-16.589 


21.187 


23.308 


23.502 




21.378 


22.657 


23.174 


22.955 


22.996 


23.413 




AM 11 


20.290 


-122.984 




23.441 


23.689 




21.477 


22.890 


23.206 


23.162 


23.140 






AM 12 


-9.226 


1.132 




23.407 


23.867 


20.985 






23.256 










AM 13 


-8.234 


-14.431 


21.036 


22.554 


22.269 


20.734 


21.103 


22.515 


22.685 


22.594 


22.694 


22.165 


looks BSS in WFC, optical WFPC2 


AM 14 


-11.461 


2.148 




23.679 


23.867 


21.198 






23.505 










AM 15 


-2.753 


2.081 


21.066 


23.448 


23.444 


20.644 


21.226 


22.892 


24.173 


22.822 


22.885 


22.662 


almost blended with BSS in WFC 


AM 16 


3.907 


-38.774 


20.348 


21.752 


21.853 




20.301 


21.457 


21.686 


21.640 


21.704 


21.822 




AM 17 


25.687 


-57.545 


20.390 


22.599 


22.796 








22.469 










AM 18 


-9.150 


44.548 


20.530 


24.649 


24.042 


















AM 19 


-3.935 


-13.644 


20.971 


23.103 


23.193 


20.628 


21.059 


23.047 


22.996 


23.016 


22.948 


23.261 




AM 20 


3.324 


6.803 


19.497 


19.224 


19.077 








19.305 











Note. — Table is available (along with photometric errors) in its entirety via the link to the machine-readable version above. PA: Post-AGB; SH: Supra-HB; AM: AGB 
manque; V: RR Lyrac variable; A: AGB; H: HB; R: RGB. The table only includes stars from the union of the ACS WFC and WFPC2 fields for proposals 7628 and 7630. 



